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A common method of adjusting the metal-insulator transition temperature of Mi VO 2 is via 
disruption of the Peierls pairing by doping, or inputting stress or strain. However, since adding 
even small amounts of dopants will change the band structure, it is unclear how doped VO 2 retains 
its insulating character observed in experiments. While strong correlations may be responsible 
for maintaining a gap, theoretical evidence for this has been very difficult to obtain due to the 
complexity of the many-body problem involved. In this work we use GW calculations modified to 
include strong local k-space interactions to investigate the changes in band structure from tungsten 
doping. We find that the combination of carrier doping and the experimentally observed structural 
defects introduced by inclusion of tungsten are consistent with a change from band insulating to 
Mott insulating behaviour. 

PACS numbers: 71.30.+h,71.27.+a,74.20.Pq,75.10.-b,71.20.-b 


INTRODUCTION 

Strongly correlated materials are an enormously ac¬ 
tive area of research as these systems exhibit many useful 
phenomena [I], with high temperature superconductivity 
being the most obvious example. However research into 
such compounds has suffered in comparison to devices 
based on weakly correlated materials due to the inappli¬ 
cability of the rigid band model of electronic structure to 
such compounds For weakly correlated materials, 

in which electron-electron interactions are a very small 
perturbation, band theoretical techniques such as Den¬ 
sity Functional Theory (DFT) have proven to be accurate 
enough predictors of material behaviour to significantly 
facilitate device design 3], In strongly correlated mate¬ 
rials however, band theory fails spectacularly [2], such 
that even qualitative predictions are almost always in¬ 
accurate. For example, DFT characteristically predicts 
Mott insulators to be metallic 0®, a significant issue 
when designing materials to take advantage of their in¬ 
sulating behaviour. 

Vanadium dioxide is considered a prototypical strongly 
correlated material, which undergoes a first order in¬ 
sulator metal transition at ~ 340 K from a low tem¬ 
perature monoclinic structure to a high temperature, 
metallic tetragonal structure [5|. In the last decade, ad¬ 
vances in nanofabrication techniques have seen an explo¬ 
sion of interest in devices based on VO 2 nanobeams mm, 
which show enormous potential for sensors [8] and even 
new transistor gates [3]. Dynamical Mean Field The¬ 
ory L 10j and GW calculations m suggest that pure VO 2 
is a band insulator resulting from the three-dimensional 
Peierls pairings occurring in the transition from tetrago¬ 


nal to monoclinic, with strong correlations in the metallic 
form responsible for driving the transition to the insulat¬ 
ing state m- 

However, many methods of controlling the critical tem¬ 
perature revolve around disrupting this Peierls pairing, 
such as the input of stress or strain [12], or doping [15] . 
Doping with tungsten has long been known to reduce 
T c by ~ 23 K per atomic percent [TJ. However, when 
viewed through the lens of band theory, this presents 
something of a paradox, as doping carriers into pure VO 2 
should result in a metallic structure. Experiments con¬ 
firm however, VO 2 doped with < 10 % tungsten remains 
insulating M ■ In addition, photoemission experiments 
reveal that tungsten assumes the place of a vanadium 
atom in the VO 2 lattice, but does not Peierls pair with 
neighbouring atoms |14| , despite having the valence elec¬ 
trons to do so (unlike e.g. doping with Ti). Rather, the 
local environment of the tungsten is tetragonal. 

Figure 1 compares the different geometries involved: 
Figure la illustrates the Peierls paired and antiferroelec- 
trically distorted vanadium chain in the undoped com¬ 
pound. Figure lb illustrates a substitutionally doped 
chain, in which the tungsten adopts the position of a 
vanadium atom. Figure lc corresponds to the exper¬ 
imentally observed |141 Il5j tungsten environment; the 
tungsten atom sits equidistant from the neighbouring 
vanadium atoms, and in the center of the oxygen oc- 
tadehron, i.e. no Peierls pairing or antiferroelectric dis¬ 
tortion. Given that the insulating state occurs via the 
vanadium atoms in the tetragonal structure forming pairs 
[5], disruption of this structural rearrangement as per 
Figure lc, combined with the extra carriers would ap¬ 
pear from a band theoretical view to transform tungsten- 
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doped VO 2 to a metallic structure. Unless of course 
the experimentally observed insulating behaviour in the 
doped structure arises due to strong correlations in the 
partially filled band. 

The resolution of such a paradox from a computational 
perspective has been out of reach, due to the failure of 
ab initio methods to reproduce Mott physics [3]. Re¬ 
cently however, it was revealed that GW calculations can 
be modified to include strong local k-space correlations 
(the “Partially Screened GW” method or PS-GW US]), 
and can reproduce two significant characteristics of Mott 
physics: the splitting of partially filled bands into the 
upper and lower Hubbard bands, and the giant transfer 
of spectral weight upon carrier doping. The ability of 
this approach to reproduce these characteristics means 
that it is now possible to investigate systems in which 
the rigid band model fails, and in this work we apply 
it to the study of tungsten-doped VO 2 . We find that 
the non-stoichiometric number of electrons results in a 
transition from band to Mott insulating behaviour. 

METHODS 

The Mi VO 2 structure El was first relaxed to the 
ground state using PBE-GGA [T5] DFT jT5] and the 
Brillouin zone integration method of Bloechl et al. j20] ■ 
6 x 6 x 6 , 4x6x6, 4x4x4 Monkhorst-Pack [21] k-space 
grids were used for the 25 %, 12.5 % and 3.1 % struc¬ 
tures, which were comprised of one, two and eight VO 2 
unit cells in lxlxl, 2 x 1 x 1 , and 2 x 2 x 2 configura¬ 
tions respectively, with a single vanadium atom replaced 
with tungsten. The GW calculations were performed us¬ 
ing the implementation of Shishkin and Kresse [22] [233 
as contained in the Vienna Ab Initio Simulation Package 
(VASP) (23|, after first calculating input wavefunctions 
using DFT with PBEjJS] functionals. The standard GW 
calculations were performed on a frequency grid of 30 
points, using a cutoff energy of 200 eV. Strong correla¬ 
tions were included by setting the derivatives of the wave- 
functions with respect to the k-point grid to zero (the PS- 
GW method, see Booth et al. mi). Five self-consistency 
steps were used for the 12.5 % doped structures, while 
four and two steps were used for the 12.5 % and 3.1 % 
doped structures respectively, and all strongly correlated 
calculations were performed at a single frequency point, 
u >=0 m • In all calculations the Fermi level is set to zero 
energy. 

RESULTS AND DISCUSSION 

We investigate this problem by first examining the ef¬ 
fect of substitutionally exchanging one vanadium atom 
for tungsten on the band structure and density of states 
(DOS) of VO 2 , and then reduce the concentration to 


more experimentally accessible values. Figure 2a illus¬ 
trates the results of unmodified G 0 W 0 calculations on 
the Mi structure of VO 2 (corresponding to the vana¬ 
dium chain environment of Figure la). The GW band 
structure (black filled circles) clearly exhibits splitting of 
the valence and conduction bands at E^ in contrast to 
the DFT bands (gray lines), which do split, however the 
gap is almost non-existent. The calculated GW gap is ~ 
0.70 eV, in excellent agreement with photoemission (also 
~ 0.70 eV)|25|. Therefore, for Mj VO 2 , screening the 
Hartree-Fock interaction with the dielectric matrix re¬ 
sulting from non-interacting Green functions is sufficient 
to reproduce the experimental characteristics of the un¬ 
doped structure. 

However, Figure 2b indicates that for substitutionaly 
tungsten-doped VO 2 (Figure lb) this is not the case. The 
red trace corresponds to the DOS of an unmodified G 0 W 0 
calculation of a VO 2 unit cell with one tungsten atom 
replacing one of the vanadiums. The DOS suggests that 
upon doping, the Fermi level shifts into the conduction 
band, which is dragged down slightly, reducing the gap 
between the original valence band and conduction band. 
The shift of the Fermi level results in a metallic structure, 
with a large number of empty states sitting at and just 
above Ep’. This is directly contradicted by experimental 
data suggesting that tungsten-doped VO 2 is insulating 
in the monoclinic form M- 

This situation embodies the failure of the rigid band 
model for strongly correlated materials. Pure VO 2 has 
a filled valence band (Figure 2a). Therefore, any doped 
states will necessarily inhabit the conduction band. How¬ 
ever the conduction states consist of a broad, continuous 
mix of bands. There are no gaps, and therefore shifting 
the Fermi level into it will create a metallic structure. 
Figure 2c illustrates the band structure corresponding 
to the G 0 W 0 data of Figure 2 b. Obviously, the “rigid” 
bands result in the Fermi level intersecting the conduc¬ 
tion band minima at T, C, E and D. The good agreement 
between the G 0 W 0 and DFT data of Figure 2c makes it 
apparent that there will be a continuum of states avail¬ 
able above Ej?, as the bands vary smoothly in k-space. 
Therefore, since pure VO 2 has a filled valence band, ac¬ 
cording to band theory, any amount of dopant, small or 
large, will result in a metallic structure due to the con¬ 
tinuous nature of the states above E^. 

Given that tungsten-doped VO 2 is an insulator, and 
taking into account the existence of the Mott insulating 
M 2 structure of VO 2 which can be brought about by Cr 
doping 26], the most obvious place to look for an ex¬ 
planation is if the structure becomes a Mott insulator 
upon doping with tungsten. We can gain some idea of 
the plausibility of this by first transforming the doped 
structure into a Mott insulator, and then looking at the 
effects of the experimentally observed structural distor¬ 
tion in this context. Adding strong local k-space correla¬ 
tions generates the black trace of Figure 2b (PS-G 5 W 5 ), 
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and the band structure of Figure 2d. Comparing the 
DOSs of Figure 2b reveals that strong correlations split 
the filled states below from the bulk of the conduc¬ 
tion states, generating a broad low density region strad¬ 
dling the Fermi level. Comparing the band structures of 
Figures 2c-d indicates that this feature results from the 
filled conduction band minima at and around F, C, E 
and D shifting to lower energy, while the empty states 
are pushed upwards. 

Figure 3 reveals the effect of increasing the dopant 
level. The DOS data of Figure 3a shows that adding 
another tungsten atom (black line) causes more states to 
cross the gap from the conduction band, adding to the 
small peak observed in the one atom doped structure (red 
line) at ~ -1.0 eV, revealing it to be the lower Hubbard 
band. Such an effect is consistent with the transfer of 
spectral weight expected in a Mott material (2T], and ob¬ 
served in photoemission experiments on TiOCl |28j . and 
in PS-GW data of CuO [lB]. The band structure of Fig¬ 
ure 3b indicates that these states again come from the 
conduction band minima, with more states from around 
the minima at T, C, E and D crossing the gap, and 
also one of the flat bands in the A E region being 
filled. The magnitude of the DOS at Ep is roughly un¬ 
changed by doping, and considerably smaller than the 
GoWo structure (Figure 2c), suggesting that strong local 
correlations will result in a structure which is less metal¬ 
lic. However, experimentally, tungsten-doped VO 2 also 
exhibits a symmetrisation of the internuclear distances 
around the tungsten atom (Figure lc). Figure 4 investi¬ 
gates the effect of this by comparing standard GoWo data 
with a more strongly correlated PS-G 5 W 5 calculation. 

Figure 4a plots the DOS data of each calculation, and 
reveals that according to standard G 0 W 0 calculations 
(red trace) the experimentally observed symmetrisation 
of the tungsten atom results in the original band gap 
closing completely, generating a metallic structure with 
an ~ 3 eV peak in the DOS straddling Ejr. Thus, band 
theory suggests that rather than the structure rearrang¬ 
ing to preserve an insulating state, it becomes even more 
metallic. In context this is unsurprising given that VO 2 is 
metallic above T c via the adoption of a tetragonal struc¬ 
ture. Therefore, according to band theory, distorting the 
structure to reduce Peierls pairing takes the monoclinic 
form towards the tetragonal form, closing the gap. In¬ 
cluding strong correlations in the (GsWsTo) calculation 
however reveals very different behaviour. The original va¬ 
lence band and the doping-induced lower Hubbard band, 
corresponding to the band minima at T, C, E and D, 
have been combined into a single valence band, split from 
the upper Hubbard band by a “gap” of ~ 1.2 eV. 

This symmetrisation of the local environment of the 
tungsten dopant observed in photoemission experiments 
PIH5| (see Figure lc) can be reconciled by starting from 
the Mott insulating ansatz. From Figures 2b-d, the most 
significant effect of introducing a tungsten dopant is to 


dope an electron into the conduction band. From figures 
2a and 3a, it is apparent that undoped Mi VO 2 already 
has a filled valence band, therefore doping introduces an 
unpaired electron into the system, creating a partially 
filled band. Substitutional doping of a tungsten atom 
into the crystal structure (as observed in photoemission) 
will result in an atom in one of the vanadium chains hav¬ 
ing an extra electron, therefore doping creates a partially 
filled band corresponding to an extra (unpaired) elec¬ 
tron. If the tungsten atom Peierls pairs with a vanadium 
atom, the extra electron will interact with both electrons 
of the pairing, as the pairing will decrease the hopping 
energy by reduced the internuclear spacing, and creat¬ 
ing potential overlap between the atoms. However if the 
structure symmetrises, the two valence electrons on the 
tungsten atom localise on it, and the single d-electron 
of the vanadium localises on the vanadium site, reduc¬ 
ing correlations. In this respect, this symmetrisation is 
a classic Mott transition, increasing internuclear spacing 
giving rise to smaller correlations. The fact that this 
symmetrisation coincides with the environment observed 
in the high temperature tetragonal structure supports 
DMFT data which indicates that the metallic tetragonal 
structure is highly correlated [IB] . 

However, the W 0 . 25 V 0 . 75 O 2 doped structure of Figure 
4c corresponds to a dopant percentage of 25 %, which 
is much higher than those commonly used to adjust T c . 
The data of Figure 4a suggests that there still exist states 
interstitial to the Hubbard bands, and thus this struc¬ 
ture may be a “bad metal” rather than an insulator [59]. 
Figure 4d explores the effect of reducing the dopant con¬ 
centration from 25- through 12.5- to 3.1 % (scaled such 
that they integrate to the same total number of states). 
As the dopant percentage decreases, the density of states 
at Ep decreases concurrently, suggesting that the struc¬ 
ture becomes increasingly more insulating as the dopant 
amount tends to zero, consistent with the fact that fewer 
carriers are being doped into the structure. This data in¬ 
dicates that at dopant percentages designed to drop T c 
to approximately room temperature (~ 2 %), tungsten 
doped VO 2 will exhibit insulating character. 

CONCLUSIONS 

In summary, by including strong correlations in GW 
calculations on tungsten-doped VO 2 it becomes evident 
that doping carriers into the band insulating VO 2 struc¬ 
ture results in a switch to Mott insulating behaviour. 
The added electron inhabits the conduction band min¬ 
ima, and strong local correlations in k-space split these 
minima from the rest of the bands, while a concurrent dis¬ 
ruption of the Peierls pairing created by the dopant atom 
results in them combining with the original valence band 
into a broad, mixed Peierls/Hubbard band. The data 
suggests that any modification of the pure VO 2 structure 
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which results in an increase in electron number will create 
a Mott insulating defect band in which the Peierls pairing 
is disrupted, and the structure seeks to minimise correla¬ 
tions by symmetrising the hopping distances around the 
defect. Therefore, non-stoichiometric structures of VO 2 
will exhibit a mixed band/Mott insulating character. 
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a Experimental Vanadium Chain 



b Substitutional Doping 
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c Experimental Doped Configuration 



FIG. 1: a) Undoped vanadium chain running along the 
monoclinic a-axis, b) substitutionally doped chain used 
in the calculations of Figures 2 and 3, and c) 
experimentally determined doped geometry used in the 
calculations of Figure 4. 
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FIG. 2 : a) Left panel: Go Wo band structure of Mi VO2 (black circles) and corresponding DFT bands (gray lines), 
Right panel: corresponding G0W0 density of states (DOS) b) DOSs of VO2 doped with one tungsten atom using 
G 0 W 0 (red) and PS-G5W5 (black), c) Left panel: G 0 W 0 (black circles) and DFT (gray lines) band structures of 
tungsten doped VO2, Right panel: corresponding G0W0 DOS, and d) Left panel: PS-G5W5 (black circles) and DFT 
(gray lines) band structures of tungsten doped VO2, and Right panel: corresponding PS-G5W5 DOS. 
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FIG. 3: a) PS-G 5 W 5 DOSs of VO 2 doped with one (red) and two (black) tungsten atoms, b) Left panel: PS-G 5 W 5 
(black circles) and DFT (gray lines) band structures of VO 2 doped with two tungsten atoms and Right panel: 
corresponding PS-G 5 W 5 DOS. 
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FIG. 4: a) GoWo (red) and PS-G 5 W 5 (black) DOSs of tungsten doped VO 2 with internuclear symmetrisation 
around the dopant atom, b) Left panel: GoWo (black circles) and DFT (gray lines) band structures of the 
symmetrised, doped structure, Right panel: corresponding GoWo DOS, and c) Left panel: PS-G 5 W 5 (black circles) 
and DFT (gray lines) band structures of this structure, Right panel: corresponding PS-G 5 W 5 DOS, and d) 
comparison of the DOSs of V0 2 with dopant percentages of 25 % (red), 12.5 % (black) and 3.1 % (filled blue curve). 























